m 
o 
o 

(N 



o 
o 



Condensation and Evaporation of Mutually 
Repelling Particles : Steady states and limit 

cycles 



Tapati Dutta 1 , Nikolai Lebovka 2 and S. Tarafdar 3 
1 Physics Department, St. Xavier's College, 
Kolkata 700016, India 
2 F. D. Ovcharenko Biocolloid Chemistry Institute 
42, Vernadsky Av., Kyiv, Ukraine 
3 Condensed Matter Physics Research Centre, 
Physics Department, Jadavpur University, 
Kolkata 700032, India. 



g ; Email: sujataOjuphys.ernet.in (S Tarafdar) 

Abstract 

\ We study condensation and evaporation of particles which repel each other, using a simple 
set of rules on a square lattice. Different results are obtained for a mobile and an immobile 
^ ■ surface layer. A two point limit cycle is observed for high temperature and low pressure in both 

. cases. Here the coverage oscillates between a high and a low value without ever reaching a 
^ \ steady state. The results for the immobile case depend in addition on the initial coverage. 

\ Keywords: Condensation and evaporation, computer simulation, limit cycle 
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1 Introduction 

O ! Recently computer simulations on cellular automata type models [1, 2, 3]have been suggested 
| which can be applied to simulate the development of a surface layer under different physical 
# ^ ■ conditions. A randomly occupied square lattice is taken and particles are allowed to stick at 
^ ■ vacant sites or evaporate from occupied sites according to a prescribed set of rules. Appropriate 
probabilities for the 'growing' and 'culling' are set according to temperature and pressure 
conditions required. We can also introduce an effective interaction between the particles by 
making the probabilities a function of the surroundings of a site, e.g. the number of nearest 
neighbours. 

In our earlier work [1] we started with a randomly occupied square lattice and let particles 
evaporate with a probability which decreased with the number of occupied neighbouring sites. 
This amounted to an attractive interaction between the particles. The condensation was a 
function of only the superincumbent pressure. The present paper is an extension of the earlier 
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work, here we introduce an effective repulsive interaction between particles, by allowing the 
condensation probability at an empty site to decrease with the number of occupied neighbours. 
The evaporation probability is a function of temperature only. The results are very interesting. 
We study as in the previous work [1] evolution of a mobile layer by analytical and numerical 
methods and an immobile layer by computer simulation. In both cases we find a steady final 
coverage under some conditions, but a two-point limit cycle showing oscillations in the coverage 
when temperature is high and pressure low. 

In the next section we describe the model. The results are presented in the third section 
and discussion and conclusions in the last. 



2 The Model 

The process of condensation/evaporation or adsorption/desorption is described traditionally 
by two different sets of models - one for a mobile adsorption layer and one for an immo- 
bile layer [4, 5]. Probabilities for sticking and evaporation on a two-dimensional monolayer are 
specified according to the physics behind the model. These are functions of the temperature, su- 
perincumbent pressure and the existing coverage. At equilibrium, the sticking and evaporation 
probabilities are set equal and the resulting equation is solved to get the equilibrium coverage 
at that temperature and pressure. 

The simplest mobile layer model is a two-dimensional ideal gas, and the improved versions 
include interaction between particles, similar to a two-dimensional Van der Waal gas. The so- 
called immobile layer models introduce a sticking probability, depending on how long a molecule 
in the vapor above the surface is in contact with a surface site. The simplest 'immobile model' 
is the Langmuir equation derived as follows. 

For vapor condensation the particle flux, i.e. the number of particles deposited per unit 
time per unit surface area is equal to: 

where P is the pressure, 



A = ^h 2 /{2TimkT) 

is de Broglie length, h is Planck constant. The probability of condensation per unit time is 
obtained by dividing c by the total number of sites available for condensation. And the same 
may be done for the evaporation probability. Flux of evaporation from saturated surface (at 
p — 1) may be approximated by [6] : 

kT 

d=—eM-E e /kT) (2) 
Condensation probability is set equal to evaporation probability, giving equilibrium: 
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PA 3 

(l-p)=exp(-E e /kT)p (3) 



kT 

and we have the simple Langmuir equation 

pb= T^J, < 4 > 

where 

A 3 

b = — exp(P e //cT), 

E e is activation energy for evaporation. 

In this approximation it is assumed that evaporation energy is the same for all configura- 
tions. 

In the present work we introduce a repulsive interaction between particles, as follows - 
condensing particles are less likely to stick at sites surrounded by more neighbors. This mimics 
a situation where the adsorbing particle has to overcome the potential barrier due the repelling 
neighbors to reach a vacant site on the substrate. It is implied that the substrate attracts the 
adsorbing particle. 

We define a probability for condensation on a vacant site f c determined by the pressure as 

fc = P/Ps (5) 

where P s is a saturation pressure and < f c < 1. The growth probability at a vacant site with 
n occupied neighbors is proportional to f c n . We have 

Pgr = (1 - P?f* + 4(1 - p) 4 pf c 3 + 6(1 - p) Yfc + 4(1 - p)Vfc + (1 - P)P 4 (6) 

The successive terms in the RHS above represent the probability of a vacant site having re- 
spectively 4,3,2,1 and vacant near neighbor sites. P gr can be written in the compact form 

P gr =(l-p)(pf c + l-p) 4 (7) 

The culling probability is assumed to be a function of temperature only 

P cull = P (kT/P s X 3 )exp(-E e /kT) (8) 

P s is temperature dependent, and this must be taken into account. But we can approximate 
it by the simple relation that follows from equilibrium relation on the surface of fluid 

kT 

P s = —exp(-E vap /kT) (9) 
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where E vap is the enthalpy of vaporization of a fluid that supports vapour. 

In this formulation the simplest condensation evaporation equilibrium condition is 

P kT 

p(l-ri = ^ ex P (-E e /kT)p = exp -AE/kTp (10) 
where AE = E e — E vap 

Here too condensation probability f c (or reduced pressure P*) is 

fc = P/Ps (11) 

and the evaporation probability is 

U = exp(-AE/kT) (12) 

That means that the molecules evaporate from the our surface with probability 1 when 
E e = E vap , when energy of evaporation from surface E e equals to energy of evaporation from 
surface of liquid that support vapor E vap . 

So, we can rewrite P cu u as 



P cull =pexp(-1/T*) (13) 

where f c = P/P s and T* = kT / AE 

Our model is complementary to the Fowler-Guggenheim model [4, 5], where the condensa- 
tion probability is dependent only on pressure, but the evaporation probability decreases with 
the number of occupied near neighbors. The neighbors are assumed to attract the particle and 
create a barrier to evaporation. This type of interaction was incorporated in Dutta et al [1]. 

Now we study the behavior of this system for different values of the temperature and pressure 
parameters. We look at the two different situations the mobile surface layer and the immobile 
layer. 



2.1 Mobile surface layer 

Let us suppose that the energy barrier for motion of the deposited particles along the surface 
is very low so a continuous rearrangement is going on. In this case the distribution of particles 
always becomes random and the probabilities in eq(6) are always valid. So one can equate the 
growth and culling probabilities in the steady state and solve for the coverage p. We confine 
ourselves to ranges of T* and P*, which give growth and culling probabilities in the range - 
1. We find that meaningful positive values for the coverage between and 1 are obtained for a 
very narrow range of temperature and pressure. The results are shown in the next section. We 
may also see the how the coverage evolves using the following iterative procedure 

Pf=Pi + Pgr - Pel (14) 
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and 

Pi = Pf (15) 

here pi is the initial coverage and pf the final. If there is a steady solution the final coverage 
Pf converges to a steady value. We find however that in some combinations of the pressure 
and temperature parameters the coverage converges to a stable limit cycle, oscillating between 
two positive and realistic values. The solution for p in eq.(6) now becomes an unstable fixed 
point. This is very interesting and reminiscent of the discovery of oscillatory behaviour in real 
chemical systems [7, 8, 9]. 

The results obtained here can be verified from a computer simulation on a finite system. 
This was done as follows. 

In the mobile case, the adsorbed molecules can move laterally on the surface. The physical 
situation simulated is described by eq.(6) . A two-dimensional square lattice of unit spacing 
and size 300 x 300, is occupied randomly with an initial coverage Pinitiai- Every occupied site is 
assigned the value 1, empty sites are assigned the value 0. The occupied sites are then culled 
parallally with a probability which is the same for each occupied site 

g = exp(-l/T*) 

. For very large T* the culling probability approaches 1. The vacant sites are filled with a 
probability dependent on the number of vacant near neighbor sites n as fj l . After one round 
of growth and culling is complete, the concentration of the occupied sites p final is calculated. 

In the next time-step, the Pfi na i of the previous time-step becomes the new Pinitiai- The 
square lattice is then randomly occupied afresh with this Pinitiai- A complete time-step begins 
with the random occupation of all sites with a Pinitiai and ends with the assignment of the Pfi na i 
to the Pinitiai of the next time-step. This iterative process stops when pfi na i saturates with 
increasing time to a definite value. In the simulation, we checked up to 50,000 time-steps. The 
'mobility' of the molecules is simulated by the randomization of the concentration Pmitiai in the 
beginning of every time-step. 

2.2 Immobile Surface Layer 

If we consider the energy barrier to surface movement of the particles to be large enough the 
arguments presented in the last subsection are not valid. We start with a certain concentration 
of particles randomly distributed on a square lattice and allow growing and culling according 
to our rules. But after one time-step the distribution is no longer random, so we cannot find 
out the number of new growth or evaporation sites as in eq.(6). We can however simulate the 
system as above without the randomization procedure after each round of growth and culling. 
This mimics the development of a surface layer where the particles are not free to move laterally 
on the surface. The results in this case are quite different . 
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3 Results 



We present the results for the mobile and immobile layer case. The mobile layer results are 
done analytically or by numerical iteration and verified by computer simulation. The immobile 
layer results are done by computer simulation as calculation is not possible here. 

3.1 Results for the mobile layer 

The results are shown in figure 1. The final coverage p for each P*, on increasing T*, from a 
low value shows a bifurcation at a certain point to a two state limit cycle. The temperature 
and coverage at this point are denoted respectively by T*t> and pb, the inset shows the variation 
of T* b and p b with P*. 

These results are independent of the initial coverage. For T* very large and P=l there is 
an oscillation between Pimtiai and 1 — Pmitiai, as is also evident from eq.(6). 

3.2 Results for the immobile surface layer 

The immobile situation can be studied only by computer simulation. Starting with an initial 
random distribution of particles the system is allowed to evolve and the final coverage for 
different P and T* is studied. In this case as there is no randomization, the system takes a very 
long time to reach the steady state, and we have observed up to 30,000 time steps. Steady state 
solutions of p final are obtained for P cuU < 1. The steady state isobars and isotherms are shown 
in Fig. 2. and Fig.3. respectively. At P cu u=1.0, oscillations between a two-point limit sets in 
for all values of P. However for every P, there is a certain initial coverage, p a that yields a 
final coverage oscillating in a narrower confine between 0.4 and 0.5 for all P. Even a departure 
of 0.05 from p a , shifts the final coverage to the fixed point limits of the wider limit cycle. The 
rate of convergence to p final from different Pinitiais other than p a , is slower for higher values of 
P. For pmitiai very close to p a in the case of P = 0.9, we checked for 50,000 time steps to detect 
whether the final coverage was approaching the fixed points of the wider limit cycle. The p a 
increases slowly with P as displayed in Tabled. 

4 Discussion 

We have demonstrated that the condensation/evaporation problem seen as a modification of 
the percolation problem on a square lattice shows very interesting possibilities. There is firstly, 
the question of mobile and immobile layer models of chemical physics [4]. From our point of 
view equating the growth and culling probabilities to find the steady state coverage is valid 
only when the particles can rearrange very fast to preserve the random distribution. And we 
demonstrate that making the particles truly immobile,on the two-dimensional surface gives a 
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quite different final coverage. This is also dependent on the initial coverage. This was discussed 
in the previous work [1] with an attractive interaction between particles. Careful experiments 
are needed to test the justification of these arguments. 

The most striking result of this work is the observation of oscillations in the condensa- 
tion/evaporation problem, using such an extremely simple model, without taking any details 
of the system into account. Earlier well known work on non-linear mathematical equations [7] 
show the appearance of bifurcations, repeated period doubling sometimes leading to chaos. And 
such situations are found in real life in problems of physics, chemistry and biology. Chemical 
systems are known to exhibit oscillations in concentration of reactants and products in chemical 
reactions [8, 9]. There is also some study of such a phenomenon in condensation/evaporation 
problems [10, 11]. Experiments with repelling particles adsorbing on a substrate are reported 
[12], but in this experiment two different types of particles are involved. 

One possibility of having a repulsive interaction is with similarly charged particles condens- 
ing on an oppositely charged substrate. 

So in conclusion the simple stochastic model with growth and evaporation on a square 
lattice gives rise to a host of interesting phenomena. With an attractive interaction between 
particles phase transitions and hysteresis was observed [1], While a repulsive interaction gives 
oscillations in coverage under certain conditions. 

References 

[1] T Dutta, N I Lebovka and S Tarafdar, Phys Rev E 66, (2002), 046104 

[2] S S Manna, T Datta, R Karmakar, S Tarafdar, IJMPC, 13, (2002), 319. 

[3] J A GurneyE A Rietman, M A Marcus, M P Andrews, J Alloys Comp, 290, (1999) 216. 

[4] Surface Chemistry, E M McCash {Oxford University Press, New York 2001) 

[5] On Physical Adsorption, S Ross and J P Olivier {Inters cience, John Wiley and Sons, New 
York 1964) 

[6] Chemistry of interfaces, M J Jaycock and G D Parfitt, {Wiley, Chichester, 1981). 

[7] Chaos Hao Bai-Lin (World Scientific, Singapore 1984). 

[8] C G Takoudis, L D Schmidt and R Aris, Surf Sci 105, (1981), 325. 

[9] Z Hou, L Yang, H Xin, J Chem Phys 111, (1999), 1592. 

[10] K I Patrylak, O M Tarnookha, Zeolites 18, (1997), 7. 



7 



[11] V P Zhdanov, B Kasemo, Phys Rev E 61, (2000), R2184 

[12] B Riedmuller, I M Ciobica, D C Papageorgopoulos, F Frechard, B Berenbak, A W Kleyn, 
R A van Santen, J Chem Phys, 115, (2001),5244. 



p 


0.1 


0.2 


0.3 


0.5 


0.6 


0.7 


0.8 


0.9 


Pa 


0.23 


0.25 


0.265 


0.305 


0.33 


0.35 


0.39 


0.43 



Table I : The variation of p a with P 
Figure captions: 

Figure 1 : Isobars for the final coverage for the mobile layer. At T*^, the steady solution for p^ 
shows the bifurcation to a limit cycle. The inset shows variation of T*b and pb with P*. 

Figure 2 : Isobars for the immobile case showing variation of coverage with exp(—l/T*). The 
values of P for each curve are shown in the figure. 

Figure 3 : Isotherms for the immobile case showing variation of coverage with pressure P. 
Stable p are obtained for exp(— 1/T*) upto 0.9, for exp(— 1/T*)=1.0 oscillations are obtained. 
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